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Abstract: The community assembly rules and species coexistence have always been interested by 
ecologists. The community phylogenetic structure is the consequence of the interaction process between 
the organisms and the abiotic environment and has been used to explain the relative impact of abiotic and 
biotic factors on species co-existence. In recent years, grassland degradation and biodiversity loss have 
become increasingly severe on the Mongolian Plateau, while the drivers for these changes are not clearly 
explored, especially whether climate change is a main factor is debated in academia. In this study, we 
examined the phylogenetic structure of grassland communities along five transects of climate aridity on 
the Mongolian Plateau, and analyzed their relations with environmental factors, with the aims to 
understand the formation mechanism of the grassland communities and the role of climatic factors. We 
surveyed grassland communities at 81 sites along the five transects, and calculated their net relatedness 
index (NRI) at two different quadrat scales (small scale of 1 m? and large scale of 5 m°) to characterize the 
community phylogenetic structure and analyze its relationship with the key 11 environmental factors. We 
also calculated the generalized UniFrac distance (GUniFrac) among the grassland communities to quantify 
the influence of spatial distance and environmental distance on the phylogenetic 8 diversity. The results 
indicated that plant community survey using the large scale quadrat contained sufficient species to 
represent community compositions. The community phylogenetic structure of grasslands was significantly 
overdispersed at both the small and large scales, and the degree of overdispersion was greater at the large 
scale than at the small scale, suggesting that competitive exclusion instead of habitat filtering played a 
major role in determination of community composition. Altitude was the main factor affecting the 
community phylogenetic structure, whereas climatic factors, such as precipitation and temperature, had 
limited influence. The principal component analysis of the 11 environmental factors revealed that 94.04% 
of their variation was accounted by the first four principal components. Moreover only 14.29% and 
23.26% of the variation in community phylogenetic structure were explained by the first four principal 
components at the small and large scales, respectively. Phylogenetic 8 diversity was slightly significantly 
correlated with both spatial distance and environmental distance, however, environmental distance had a 
less explanatory power than spatial distance, indicating a limited environmental effect on the community 
phylogenetic structure of grasslands on the Mongolian Plateau. In view of the limited effect of climatic 
factors on the community phylogenetic structure of grasslands, climate change may have a smaller impact 
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on grassland degradation than previously thought. 
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1 Introduction 


The formation and maintenance of community biodiversity have always been the focuses of 
community ecology research (Roughgarden, 1983; Tilman, 2004; Zhang et al., 2017). 
Conservation community ecologists tend to treat all plant species equally and ignore the 
differences among species with different origins and evolution, which may lose much important 
information (Vamosi et al., 2009; Marc et al., 2012). The development of community 
phylogenetic structure, however, provides a better understanding of the community-assembling 
mechanisms (Webb et al., 2002; Swenson et al., 2007). Generally, closely related species tend to 
share similar niches (Burns and Strauss, 2011). Therefore, if community composition is 
dominated by habitat filtering, the closely related species will be gathered, and local communities 
exhibit phylogenetic clustering. Moreover, in communities where competition greatly affects 
composition, distantly related species are expected to co-occur more frequently, and as a result, 
the community is phylogenetically overdispersed. If neutral assembly is the primary factor or 
habitat filtering works with competitive exclusion, the community may be phylogenetically 
random (Webb et al., 2002; Cavender-Bares et al., 2004). The community phylogenetic structure 
could thus be a practical indicator of understanding how biotic and abiotic factors influence the 
community assembly (Webb et al., 2002; Donoghue, 2008). 

The impact of global climate change on terrestrial ecosystem degradation has attracted much 
attention (Melillo et al., 1993; Cramer et al., 2001). However, the main reason for grassland 
degradation, either climate change or human interference, has been controversial in academia 
(Cnr et al., 2008; Chen et al., 2014; Hilker et al., 2014; Cao et al., 2015). The study of community 
phylogenetic structure could provide a potential evidence about how climate change influences 
community and species assembly in grasslands (Willis et al., 2008; Lavergne et al., 2010). For 
example, a community where abiotic filtering is the dominant force should be more sensitive to 
climate change (Willis et al., 2008). In contrast, a community where competitive exclusion is 
dominant would be expected to be more stable and have a greater capacity to adapt to the changed 
climate (Cadotte et al., 2012). Thus, the relative importance of abiotic filtering in community 
assembly could be used to identify how climate and climate change influence the local 
community (Buckley and Kingsolver, 2012). 

Previous studies on community phylogenetic structure mainly concentrated on forest 
ecosystems (Kembel and Hubbell, 2006; Letcher, 2010; Baraloto et al., 2012), whereas fewer 
studies have been reported on grassland ecosystems. Some researchers have reported the 
importance of competitive exclusion in grassland community assembly (Slingsby and Verboom, 
2006; Qian et al., 2014). For example, a study of 11 Stipa grasslands in Spain revealed that most 
regions were phylogenetically random due to the collective effect of habitat filtering and biotic 
interactions within the community (Soliveres et al., 2012). However, the interactions between 
species (especially competition) were the main driver of community assembly. Yang et al. (2018) 
found that nitrogen addition significantly changed the community structure from being 
phylogenetically random to overdispersed by changing the interspecific competition and the 
colonization of species, whereas water addition only had a little influence. Studies also showed 
that community phylogenetic structure was strongly influenced by the spatial scales (Kraft and 
Ackerly, 2010; Soliveres et al., 2012). For instance, the community phylogenetic structure of 
some tropical rainforests tends to be clustered at a small scale (i.e., <100 m°) and overdispersed at 
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a large scale (Swenson et al., 2007). However, little is known on the influence of special scale on 
grassland community phylogenetic structure. Meanwhile, research on community phylogenetic 
structure of grasslands at large spatial scales (regional to global) is lacking. 

The Mongolian Plateau, located in arid and semi-arid areas, is well known for its broad natural 
grasslands, and plays a crucial role in ecological stability and biodiversity not only in North 
China but also in the whole of Northeast Asia (Chen et al., 2007). In recent decades, the 
grasslands on the Mongolian Plateau have experienced severe degradation and biodiversity loss 
(Zhang et al., 2011), and their conservation is urgently needs (Zhao et al., 2015). 

In present study, based on the grassland survey and analysis of their community phylogenetic 
structures along gradients of climate aridity on the Mongolian Plateau, we analyzed the 
relationship between community phylogenetic structure and environmental factors with the aim to 
quantify the relative importance of abiotic filtering in community assembly. Furthermore, we also 
calculated the phylogenetic B diversity among the grassland communities and analyzed their 
relations with environmental factors and spatial distance. We hypothesized that competitive 
exclusion rather than habitat filtering dominates the community phylogenetic structure of 
grasslands on the Mongolian Plateau, and that changes in climatic factors does not significantly 
influence the community phylogenetic structure of grasslands. The research of community 
phylogenetic structure of grasslands could provide theoretical and practical guidance for the 
conservation of grassland ecosystems on the Mongolian Plateau as a unified entity. 


2 Materials and methods 


2.1 Study area 


The Mongolian Plateau is located in the hinterland of Asia, and consists of the whole territory of 
Mongolia, most parts of the Inner Mongolian Autonomous Region of China and the southern 
Russia (Zhou et al., 2017). The central of this area is comprised of high plains and its edge is 
comprised of moutains, with the Da Hinggan Mountains to the east, the Altay Mountains to the 
west, the Yinshan Mountains to the south and the North Sahyan-Hentiy Mountain line to the north 
(41.6°-52.2°N, 87.6°-119.9°E; 1000-1500 m a.s.l.; Lu et al., 2009). The Mongolian Plateau is 
characterized by a typical temperate continental climate (Bao et al., 2014). The mean annual 
precipitation on the Mongolian Plateau is generally less than 200 mm and is mainly concentrated 
in summer. However, mean annual precipitation can reach up to 400 mm in the east. The annual 
mean temperature is approximately 1.9°C (Zhou et al., 2017). The dominant vegetation type is 
steppe grassland, which can be classified into meadow steppe, typical steppe and desert steppe on 
the basis of precipitation. Furthermore, a large area of desert is distributed in the west and 
southwest (Vandandorj et al., 2015). 


2.2 Experimental design 


Four transects (named a, b, c and d; Fig. 1) were set up and surveyed in 2015 from the desert 
steppe region to the typical steppe region along the precipitation gradients. In 2016, another 
transect (named e; Fig. 1) was set up along the east-west direction, where precipitation is 
approximately 240 mm and dominant vegetation type is the typical steppe. Along each transect, 
several plots (15, 14, 8, 9 and 35 plots on the transects a, b, c, d and e, respectively; 
approximately 50 km between each closest pair) were selected in wide and flat terrain with less 
interference with a total of 81 plots. A GPS receiver was used to record the latitude, longitude and 
altitude of each plot. A 10 mx10 m square was designed in each plot, and then 5 quadrats (1 mx1 
m) were set in the square at the four vertices and the center (Fang et al., 2009). All the species and 
their individual abundances in each quadrat were recorded, and green plants were clipped and 
collected at the same time. Plant materials were oven-dried at 75°C for 24 h and then weighed. 
Each quadrat (1 mx1 m) was regarded as a representative of the local community at the small 
scale (area of 1 m°). Then, the 5 small scale quadrats were combined as one to represent the local 
community at the large scale (area of 5 m°). Additionally, 2 transects (a and b) were surveyed 
using the nested quadrats (at nine scales of 1 mx1 m, 1 mx2 m, 2 mx2 m, 2 mx4 m, 4 mx4 m, 4 
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mx8 m, 8 mx8 m, 8 mx16 m and 16 mx16 m, respectively) to ensure most of local species were 
included (Stohlgren et al., 1995). The number of species in each quadrat on the transects a and b 
at the small and large scales was compared with that recorded in the nested quadrats at the nine scales. 
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Fig. 1 Location of the surveyed plots on the five transects (a, b, c, d and e) on the Mongolian Plateau. Arrows 
along the transects a, b, c and d indicate the precipitation gradients from high to low. 


2.3 Environmental data 


Fifty-year average values of 19 meteorological factors were downloaded from WorldClim 
(http://www.worldclim.org; Hijmans et al., 2005). The spatial resolution was 30 arc seconds (ca. 1 
km). Meteorological factors with a high correlation coefficient (*>0.95) were removed and then 
10 factors were remained. Then altitude was added to this group. Finally, 11 environmental 
factors (Table 1) were selected and analyzed using multiple linear regression (MLR) models. 


Table 1 Eleven environmental factors used in this research 


Factor Description 
AMT (°C) Annual mean temperature 
MDR (°C) Mean diurnal range (mean of monthly (maximum temperature—minimum temperature)) 
TS (°C) Temperature seasonality (SD [Tavg;]) 
TAR (°C) Temperature annual range (maximum temperature of warmest month—minimum temperature of coldest month) 
IS Isothermality (MDR/TAR) 100) 
MTW (°C) Mean temperature of the warmest quarter 
MTC (°C) Mean temperature of the coldest quarter 
AP (mm) Annual precipitation 
PS Precipitation seasonality (coefficient of variation) 
PD (mm) Precipitation of the driest quarter 
ALT (m) Altitude 


Note: Tavg; indicates the average temperature of the i" month. A quarter is defined as any period of three months. SD, standard 
Bı g! p! q y p 


deviation. More details can be found at http://www.worldclim.org/bioclim. 


2.4 Phylogenetic data 


We recorded all the species in the quadrats in the online database Phylomatic 
(http://phylodiversity.net/phylomatic/) based on APG (Angiosperm Phylogeny Group) III (The 
Angiosperm Phylogeny Group, 2009). A phylogenetic supertree with branch length and 
differentiation time was regarded as outputs (Webb and Donoghue, 2005). It is worth noting that 
since the branch length of the gymnosperms on the supertree is too long and the gymnosperms 
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have a low abundance in the community (Gerhold et al., 2008; Letcher, 2010), the gymnosperms 
appearing in the individual quadrats were excluded from this study. 

Net relatedness index (NRI; Webb et al., 2002) was calculated to analyze the community 
phylogenetic structure. This index assumes that all species in all plots could form a species pool. 
First, the mean phylogenetic distances (MPDs; million years) of all species pairs in the quadrats 
were calculated. Then, the numbers of species and individual species were kept constant, and the 
species names in the species pool were selected and repeated several times randomly to obtain the 
distribution of MPDs in the random null model. Finally, the observations were standardized to 
obtain the NRI using the following formula: 


MPDos — mean ( MPDxutt) 


NRIF-1x 
SD(MPDxutt) 


(1) 


where MPDovs is the observed mean phylogenetic distance value in each quadrat (million years); 
mean(MPDnu) is the average mean phylogenetic distance value of the species obtained by 
random numbers in the phylogenetic tree (million years); and SD(MPDnpun) is the standard 
deviation of phylogenetic distance value of the species obtained by random numbers in the 
phylogenetic tree (million years). 

The NRI was calculated at both the small and large scales. Generally, if the NRI>0, the species 
in the quadrat is phylogenetically clustering. In contrast, if the NRI<0, the species in the quadrat 
is phylogenetically overdispersed. Moreover, if the NRI=0, the species in the quadrat is 
phylogenetically random, or there is no community phylogenetic structure (Webb, 2000). 

In the past, the studies of community phylogenetic structure in forest ecosystems were mostly 
based on a single species and its individual numbers in the community. However, the tussock 
graminoids prevail in the grassland ecosystems. For example, one Carex korshinskyi and one 
Stipa grandis are considered equally as one unit in the calculation of NRI, but this will obviously 
underestimate the role of S. grandis in the community. In addition, it is difficult to calculate the 
numbers of individual stoloniferous plants (i.e., Artemisia frigida). Thus, we proposed to use the 
relative dry biomass of the species in the quadrats instead of the numbers of species individuals. 
The relative dry biomass (RDB) was calculated using the following formula: 


(2) 


RDB integer Dry mass of one species 0° |. 
Total dry mass of all species 


2.5 Phylogenetic B diversity 


To further analyze the impact of environmental factors on the community composition, we 
calculated the phylogenetic B diversity. The abundance-weighted Generalized UniFrac distance 
(GUniFrac) was used to measure the phylogenetic B diversity (Chen et al., 2012). GUniFrac was 
calculated using the following equation: 


A 8 
Dr beh + pPy e | 

: dD; +D; 

GUniFrac=-——_,—____——_— > (3) 
2 bp; + p; ) 

where m is the number of branches (million years); b; is the length of the i branch (million 
years); p and p? are the branch proportions for communities A and B, respectively; and the 

parameter a is used to control the weight of abundant lineages. 


2.6 Statistical analysis 


Statistical analysis was performed using R 3.4.2 (R Core Team, 2017, Vienna, Austria). Duncan's 
test was performed using the agricolae package to analyze the difference between the species in 
the quadrats at the two scales (the small and large scales) and in the nested quadrats at the nine 
scales (1 mx1 m, 1 mx2 m, 2 mx2 m, 2 mx4 m, 4 mx4 m, 4 mx8 m, 8 mx8 m, 8 mx16 m and 16 
mx16 m) (De Mendiburu, 2014). NRI and GUniFrac were calculated using the picante package 
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(Kembel et al., 2010). The 11 environmental factors were reduced over four axes using the 
principal component analysis (PCA) to quantify the relative contributions of environmental 
factors and their interactive effects. Then the relationship between the four axes and the NRI was 
analyzed using a general linear model (GLM) to quantify the relative contribution of 
environmental factors to the NRI. Spatial distance was determined by calculating the 
geographical distance between pairs of plots using latitude and longitude. Then, the Euclidean 
distance between every two plots in the space of the first four principal component axes was 
calculated to measure the environmental distance (Qian et al., 2017). Pearson's correlation was 
used to analyze the effects of diffusion restriction (spatial distance) and habitat filtering 
(environmental distance) on phylogenetic ß diversity. 


3 Results 


3.1 Species number in the quadrats at different scales 


In our study, 331 species of angiosperms belonging to 47 families and 172 genera were recorded 
in the 81 plots. Species in the nested quadrats showed that increasing quadrat area rapidly 
increased the species numbers in the quadrat, but the increasing numbers leveled off when the 
quadrat area reached 16 m? (4 mx4 m; numbers of 22.97 (+8.72)). There was no significant 
difference of species numbers between the quadrat of 5 m? (large scale; 26.83 (+10.80)) and the 
nested quadrats of 64 m? (8 mx8 m; 27.24 (+10.23)), 128 m? (8 mx16 m; 28.31 (+10.70)) and 256 
m? (16 mx16 m; 29.48 (+11.72)) (Fig. 2). 
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Fig. 2 Species numbers in the quadrats at the two scales (S, small scale of 1 m?; L, large scale of 5 m°) and in 
the nested quadrats at the nine scales (1 mx1 m, 1 mx2 m, 2 mx2 m, 2 mx4 m, 4 mx4 m, 4 mx8 m, 8 mx8 m, 8 
mx16 m and 16 mx16 m). Different lowercase letters indicate the significant differences of species number 
among different quadrats at P<0.05 level. The upper whisker represents the maximum; the lower whisker 
represents the minimum; the top of box represents the upper quartile; the bottom of box represents the lower 
quartile; the band in the box represents the median; and the black circle represents the extreme value of species 
number among different quadrats. 


3.2 Community phylogenetic structure 


A significant community phylogenetic overdispersion was found on the Mongolian Plateau (Fig. 
3). At the small scale, there were 21 plots with positive values of NRI and 60 plots with negative 
values of NRI, and the average value of NRI was —0.20 (+0.05). At the large scale, there were 12 
plots with positive NRI values and 69 plots with negative NRI values, and the average value of 
NRI was —0.32 (+0.04). The plots with negative NRI values in the quadrats at the small and large 
scales accounted for 74.07% and 85.19% of the total 81 plots, respectively. The average NRI 
value in the quadrats at the small scale was significantly lower than that in the quadrats at the 
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large scale (P=0.049). In addition, one-tailed ¢ test indicated that the values of NRI at both the 
small and large scales were significantly less than 0.00 (P<0.01). 
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Fig. 3 Net relatedness index (NRI) values of 81 plots in grassland ecosystems on the Mongolian Plateau at the 
small scale (S, small scale of 1 m°) and large scale (L, large scale of 5 m°). Different letters indicate the 
significant differences at P<0.05 level. The upper whisker represents the maximum; the lower whisker represents 
the minimum; the top of box represents the upper quartile; the bottom of box represents the lower quartile; the 
band in the box represents the median; and the black circle represents the extreme value of species number among 
different quadrats. 


3.3 Relationship between NRI and environmental factors 


The mean annual precipitation in the study area ranged from 108 to 408 mm. However, the values 
of NRI did not change with the precipitation gradient nor with the variations of other 
environmental factors. The results of linear stepwise regression showed that the NRI at the small 
scale was significantly correlated with annual precipitation (P<0.05) and altitude (P<0.01), and 
that the NRI at the large scale was only significantly correlated with altitude (P<0.01; Table 2). 


Table 2 Multiple regression between net relatedness index (NRI) and environmental factors after stepwise 
regression 


Scale Environmental factor Estimate Standard error t P 
Small scale (Intercept) 0.3527 1.4643 0.241 0.810 
MTW 0.0337 0.0250 1.351 0.181 
AP 0.0030 0.0013 2.332 0.022* 
PS -2.0354 1.0637 -1.914 0.059 
PD —0.0440 0.0367 -1.199 0.234 
ALT 0.0005 0.0002 2.936 0.004** 
Large scale (Intercept) —1.0720 0.2326 —4.608 0.000°** 
AP 0.0011 0.0006 1.822 0.072 
ALT 0.0004 0.0001 3.121 0.003™ 


Note: MTW, mean temperature of the warmest quarter, AP, annual precipitation; PS, precipitation seasonality (coefficient of variation); 
PD, precipitation of the driest quarter; ALT, altitude. *, ** and “* indicate significance at P<0.05, P<0.01 and P<0.001 levels, 
respectively. 


Additionally, the 11 environmental factors were subjected to PCA, and the first four principal 
components accounted for 94.04% of the total variation of NRI caused by these factors. The first 
axis (PC1) was related to the temperature of the coldest quarter and the interannual temperature 
variation (i.e., temperature seasonality and temperature annual range). The second axis (PC2) was 
related to the altitude and the temperature of warmest quarter. The third axis (PC3) was correlated 
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to the precipitation of the driest quarter and the monthly temperature variation (i.e., mean diurnal 
range of temperature), and the fourth axis (PC4) was related to the annual precipitation and 
precipitation seasonality. The GLM results showed that PC3 explained the largest proportion 
(11.94% and 17.65% for the small and large scales, respectively) of the overall variation, whereas 
other axes did not make significant contributions (Table 3). All the environmental factors 
accounted for less than a quarter (14.29% and 23.26% for the small and large scales, respectively) 
of the total variation. 


Table 3 General linear model (GLM) analysis for the integrative effects of environmental factors on NRI 


Small scale Large scale 
Source 

MS SS (%) MS SS (%) 

PCl 0.02 0.16 0.33 2.92 
PC2 0.17 1.25 0.22 1.96 
PC3 1.64 11.94° 2.00 17.65" 
PC4 0.13 0.94 0.08 0.73 
Residuals 0.15 85.71 0.11 76.74 


Note: PC1, the first axis; PC2, the second axis; PC3, the third axis; PC4, the fourth axis; MS, mean squares; SS, proportion of the 
variance explained by the environmental factors; *, significance at P<0.05 level. 


3.4 Correlation of phylogenetic B diversity with spatial distance and environmental 
distance 

The phylogenetic B diversity had a weak but significant positive correlation with spatial distance 
and environmental distance (P<0.01), and the Pearson's correlation coefficients (r) between 
GUniFrac and spatial distance and between GUniFrac and environmental distance were 0.132 and 
0.050, respectively (Fig. 4). Generally, the phylogenetic B diversity (GUniFrac) increased with 
increasing spatial distance and environmental distance. However, environmental distance had less 
explanatory power for phylogenetic B diversity than that for spatial distance. 
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Fig. 4 Correlation of phylogenetic B diversity (GUniFrac) with spatial distance and environmental distance 


4 Discussion 


4.1 Community phylogenetic structure of grasslands on the Mongolian Plateau 


In this study, most of the plots had negative NRI values at both the small and the large scales 
(accounting for 74.07% and 85.19% of the total 81 plots, respectively; Fig. 3) and the mean value 
of NRI was significantly less than 0.00 at the two scales, suggesting a phylogenetic 
overdispersion on the Mongolia Plateau (Webb, 2000). In previous studies, competitive exclusion, 
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habitat filtering, and neutral assembly were used to explain the mechanisms of community and 
species assembly (Kraft et al., 2007). The significant phylogenetic overdispersion in the studied 
grasslands indicates that competitive exclusion plays a relatively important role in species 
assembly (Webb et al., 2002; Kraft et al., 2007) and that niche differentiation is the main factor in 
forming the community composition (Cavender-Bares et al., 2009) of grasslands on the 
Mongolian Plateau. 

In general, the available resources of grassland ecosystems are more limited than those of 
forest ecosystems (Chu et al., 2016). However, the main driving force of the community 
composition in this study is competitive exclusion rather than habitat filtering. 

The latter usually plays an important role in East Asian forest ecosystems (Feng et al., 2015). 
This trend may occur, because compared with woody plants, herbaceous plants tend to have shorter 
life cycles and more rapid reproduction that facilitates their occupation of suitable habitat patches 
(Donoghue, 2008). Qian et al. (2014) found that phylogenetic overdispersion dominate in the 
herbaceous ecosystem, but phylogenetic clustering is the dominate in woody assemblages in the 
same mountain ecosystem, which can be attributed to that the protection of underground rhizomes 
helps herbaceous plants resist low temperature in winter. Additionally, Ricklefs and Latham (1992) 
proposed that herbaceous plants have edaphic- and microhabitat-dominated niches rather than 
climate-dominated niches. In addition, some studies highlight a negative density dependence, which 
limits the population densities of individual species and thus has an important effect on community 
assembly at a small scale (Webb et al., 2006). The quadrat area in our study represented a relatively 
small scale, which could partially explain the community phylogenetic overdispersion. 


4.2 Effect of spatial scale on community phylogenetic structure 


Spatial scale is an important parameter that affects the structure of communities (Cavender-Bares 
et al., 2006; Swenson et al., 2006; Emerson and Gillespie, 2008). Many studies suggest that 
overdispersion is more common at small spatial scales (e.g., less than 100 m°), whereas the 
community phylogenetic structure begins to show clustering at large spatial scales greater than 
100 m’, (Kembel and Hubbell, 2006; Swenson et al., 2007). A possible reason for the clustering 
structure of communities with increasing spatial scale is the fiercer interspecific competition as a 
result of the relatively limited resources at the small spatial scale (Volkov et al., 2005). At a large 
spatial scale, the effect of habitat filtering on community composition will become more visible 
due to the greater spatial heterogeneity (Davies et al., 2005). The change in the community 
phylogenetic structure with the spatial scale is not due to the difference in the sampling scale but 
due to the environmental gradient caused by the change in the spatial scale (Willis et al., 2010). 
However, those conclusions were mainly drawn from the studies on forest ecosystems. Our 
results showed that, the NRI of grasslands tended to decrease with increasing quadrat area (Fig. 
3). Herbaceous plants generally make more investment in asexual reproduction than woody plants 
(Tamm, 1956; Hughes et al., 1988). Thus, the same species is more likely to be locally crowded at 
smaller scales (Nishitani et al., 1999). As a result, the same species is clustered around one 
mother plant, which may partially offset the phylogenetic overdispersion at the small scale. 

For this study, in view of the fact that quadrat area less than 1 m° can only contain limited plant 
species in the grassland ecosystems, these species are the embodiment of the patch in the 
community. In contrast, although the sampling area at the large scale was only 5 m°, as it is a pool 
of 5 seperated small quadrats of 1 m? within an area of 10 mx10 m, the species numbers had no 
significant difference with those in the nested quadrats ranging from 64 m° to 256 m? (Fig. 2), 
which means that the species at the large scale (5 m°) could reflect the local community 
compositions at 100 m? or larger spatial scales. Therefore, we conclude that the community 
phylogenetic structure at the large scale is a more reliable representative of the local community. 


4.3 Effect of environmental factors on community phylogenetic structure 


Temperature and precipitation are the main factors affecting the type and composition of 
communities in terrestrial ecosystems (Field et al., 1998; Fernandez, 2002). Precipitation 
especially restricts the community biodiversity and biomass on the Mongolian Plateau, which is 
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located in the arid and semi-arid area (Bai et al., 2008). In our study, although the annual mean 
temperature and mean annual precipitation vary greatly among different plots (ranging from 
—6.5°C to 6.3°C and from 108 to 408 mm, respectively), the effects of single climatic factor (e.g., 
temperature or precipitation) on community phylogenetic structure are not significant (with the 
exception of the annual precipitation having a small impact on NRI (P=0.02) at the small scale). 
Notwithstanding, our results show that community phylogenetic structure tends to be clustered at 
high altitudes (Table 2). Indeed, altitude affects the community environment, including 
hydrothermal conditions. High altitude is usually regarded as a harsher environment than low 
altitude (Hoiss et al., 2012). Therefore, the influence of habitat filtering is relatively prominent at 
high altitude, and this result is supported by Kembel and Hubbell (2006). The GLM results show 
that the third principal component of PCA (mainly related to the precipitation of driest quarter and 
mean diurnal range temperature) explains the largest proportion of the total variation in NRI. 
However, the contribution of environmental factors is still limited. Taking the first four PCA axes 
into account, the contribution of environmental factors is less than a quarter of the total variation. 
This result is consistent with the abovementioned results, indicating that habitat filtering does not 
play a major role in determining the community phylogenetic structure. Although climatic factors 
significantly influence the grassland taxonomy (Zavaleta et al., 2003; Bai et al., 2008), they may 
have different impacts on community phylogenetic structure (Webb et al., 2002). This is also 
supported by Yang et al. (2018) who revealed that water addition changes species richness but 
does not influence the community phylogenetic structure. Additionally, grasslands are influenced 
by multiple factors not only including climate change but also including human activities (Du et 
al., 2004), such as grazing, mining and urbanization, which deeply impact grassland community 
and succession (McKinney, 2002; Williams et al., 2005). Therefore, the effect of human activities 
on community phylogenetic structure warrant more attention in the future. 

Diffusion restriction and habitat filtering are considered the predominant factors in the 
formation of community B diversity (Freestone and Inouye, 2006; Kraft et al., 2011). Our study 
confirms a distance decay (in terms of both spatial distance and environmental distance) in 
grassland ecosystems on the Mongolian Plateau (Fig. 4), and the result shows that spatial distance 
has a more explanatory power for phylogenetic PB diversity than environmental distance. This 
result further suggests the limited influence of habitat filtering on community phylogenetic 
structure. However, diffusion restriction and habitat filtering only explain a small part of the 
phylogenetic B diversity, with a low correlation coefficient. This result is coincided with the 
finding of Qian et al. (2014) that environmental factors only explain a little variation in 
phylogenetic B diversity for herbaceous assemblages. This might be due to that the Mongolian 
Plateau is relatively flat with no significant geographical isolation and topographic relief present 
(Yue, 2011). In addition, the existence of widespread species also limits the variation in 
phylogenetic B diversity (Guojun et al., 2010). 

Human activity and climate change have been recognized as two main drivers of grassland 
degradation (Gang et al., 2014). Our study shows that habitat filtering does not significantly affect 
the community phylogenetic structure on the Mongolian Plateau. Therefore, the variation in 
abiotic environment along a particular gradient (i.e., increased precipitation or temperature) does 
not lead to a change in community phylogenetic structure. Furthermore, we speculate that human 
activities may play a more important role than climate change in grassland degradation. 


5 Conclusions 


Our results indicate that grasslands on the Mongolian Plateau are phylogenetically overdispersed 
and that competitive exclusion plays a more important role in the community composition of 
grasslands than habitat filtering. A negative NRI value would become more common at the large 
scale in the grassland ecosystems. Altitude has a significant positive effect on NRI, indicating a 
decline in phylogenetic overdispersion with increasing altitude. Climatic factors have a limited 
effect on community phylogenetic structure of grasslands on the Mongolian Plateau, indicating 
that climate change has a less impact on community structure of grasslands on the Mongolian 
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Plateau than previously thought. 
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